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In systems where the design Inlet and outlet pressures P am 5 are maintained 
above the thermodynamic critical pressure P c , It Is often assumed that 
two-phase flows within the system cannot occur. Designers rely on this simple 
rule of thumb to circumvent problems associated with a highly compressible 
two-phase flow occurring within the supercritical pressure system along with 
the uncertainties In rotordynamlcs, load capacity, heat transfer, fluid 
mechanics, and thermophysical property variations. The simple rule of thumb 
Is adequate In many low-power designs but Is Inadequate for high-performance 
turbomachines and linear systems, where two-phase regions can exist even 
though P,amb > Pc. Rotordynamlc-fluld-mechanlc restoring forces depend on 
momentum differences, and those for a two-phase zone can differ significantly 
from those for a single-phase zone. Using the Reynolds equation the angular 
velocity, eccentricity, geometry, and ambient conditions are varied to 
determine the point of two-phase flow Incipience, Subsequent dynamic effects 
and experimental verification remain as topics for future research. 


SYMBOLS 

a dimensionless eccentricity, e/C 
C clearance, cm 

D diameter, cm 

e eccentricity 

f,h functions scaling working fluid to reference fluid 
G r reduced mass flux, G/G* 

G mass flux, g/s-cm 2 
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' * thfermodyn^mlf parameter, (P c p c /Z c ) 1/2 , 6010 g/s cm 2 for nitrogen 

H fluid film thickness, cm 

L spacer or bearing length, cm 

l aperture (orifice, Borda) length, cm 

M molecular weight 

N rotational speed, rad/sec (rpm) 

P pressure, MPa 

P r reduced pressure, P/P c 

R radius, mm 

Rg gas constant 

S entropy, J/g K 

T temperature, K 

T r reduced temperature, T/T c 

w mass flow rate, g/s 

Z axial coordinate, cm 

Zq compressibility factor, P/pRgT 

n viscosity, Pa-s 

e circumferential coordinate, rad 

A bearing parameter 

p density, g/cm 3 

#,e shape functions 

Q angular velocity, rad/s 

Subscripts: 
amb ambient 

c thermodynamic critical point (also superscript) 

ex exit 

In Inlet 
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s saturation 

0 reference fluid 

a working fluid 


INTRODUCTION 

As a rule of thumb the extent and nature of the cavity formed downstream 
of the minimum clearance in a bearing or seal can be controlled, even elimi- 
nated, by increasing the ambient pressure. Even for systems operating above 
the thermodynamic critical pressure the rule of thumb - maintain the design 
inlet and outlet pressures well above the thermodynamic critical pressure - is 
used to circumvent two-phase flow and cavity formation (fig. 1). Uncertainties 
in stability, load capacity, heat transfer, fluid mechanics, and thermophysical 
property variations are thereby avoided. This simple rule of thumb usually 
works for low-power designs, but it is Inadequate for high-performance turbo- 
machines and linear systems. 

In many cases the two-phase flow cavity enhances bearing stability 
(refs. 1 and 2) and seal load capacity (ref. 3), The inception of cavity flow 
Is Important to both turbomachine stability and load capacity, 

Cavity flows are engendered by local "violations" of thermodynamics, 
termed "metastabllltles." Classic among these are the fine capillary work of 
Skrlpov (ref. 4) and applied work such as fracture, boiling, and condensing 
(refs. 5 and 6) . 

Thus a two-phase flow cavity within a supercritical pressure system would 
be considered unstable, perhaps even unrealistic. Such conditions do not 
violate thermodynamics but are a source of Incompatibility for the practicing 
engineer. Yet It has been demonstrated that two-phase cavity flows can be 
"nested" within a supercritical pressure system and that conditions for the 
existence of such nested cavities are prevalent for sharp-edge Inlets and for 
eccentrically loaded high-performance turbomachines (ref. 7), 

In this paper we will fix a turbomachine configuration, and for a series 
of parametric variations, establish a point of Incipient two-phase flows by 
using the Reynolds equation and methods of corresponding states. 


TURBOMACHINE CONFIGURATIONS 

Normally seals and bearings for high-performance turbomachines operate at 
small eccentricities. However, during abnormal events or at high loads, the 
dimensionless eccentricity a can approach 1 (rubbing). The load W 
supported by a stable constant-property 360° fluid film H becomes (refs. 1 
and 8) 


c -c?— 
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Departures from theory are represented by the normalized load capacity, which 
Is obtained from Integrating the numerically calculated pressure profiles 


Load capacity ratio = (W/ A AP) ca -| C /(W/A AP)^heory (4) 

Later we will make use of these relations In comparing full-film to Incipient- 
film capacities. 


Seals 


Figure 2(a) Illustrates the geometry of a high-performance turbomachine 
seal such as those used In the space shuttle main engine (SSME) (ref. 3). The 
working fluid Is parahydrogen (P c * 1.2928 MPa, T c a 32.976 K, p c = 

0.0314 g/cm 3 ). The pressure profiles (fig. 2(b)) are for a fully eccentric, 
nonrotating shaft. The Inlet pressure Is greater than P c , as Is the pres- 
sure at the Inlet to the third step. However, Immediately downstream of the 
step the pressure on one "side" of the seal (e = 0) remains above P c while 
the pressure on the opposite "side" (e ■ ir) Is nearly P s . Thus the circum- 
ferential pressure profile possesses all states from P^ n > P c at e * 0 to 
approximately P s < P c at e = ir, which represents two-phase cavity flow 
nested within a supercritical pressure regime. 

With eccentric seals pressure gradients of sufficient magnitude can exist 
to engender two-phase cavity flows both axially and circumferentially while 
embedded or nested within a supercritical pressure system. 


Bearings 

Unlike those for a seal the axial pressure gradients for a bearing are 
usually small, but the circumferential gradients are usually much larger. 
Embedded two-phase cavity flow zones In seal flows have been established and 
should occur naturally In bearings. High-speed motion pictures (to 5000 pps) 
of the Incipience, development, and collapse of the two-phase flow cavity and 
pressure-temperature profiles on the stator have been mapped for a low- 
performance machine (refs. 10 and 11). It Is also noted In reference 10 that 
the two-phase zone begins to collapse as system ambient pressure Is Increased. 
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However, for calculations using Reynolds equations (ref. 11), a high- 
performance bearing configuration In fluid hydrogen at high eccentricity 
(table I) reveals a two-phase flow cavity embedded within the system where 
p 1n > p ex > p c 3), The cavity flow may be similar to that found for 

sharp-edge Inlets (ref. 12). The pressure P drops abruptly from P^ n to 
below P s ; In the recovery zone P Increases to P s ; In the diffusion zone 
the pressure can recover or remain monotone; either case terminates In a 
two-phase choke (fig. 4). 

The onset of the embedded two-phase flow cavity zone depends on the param- 
eters of equation (1) and Is prevalent at lower ambient pressure (e.g., one- 
tenth those of table I), high eccentricity, and high rotating speed and for 
large geometries. To Illustrate the effects of two-phase flow cavity control, 
the pressure difference (Pg - P m in)/ p c 1$ tabulated (table II) as a function 
of running speed N for three eccentricities at two temperatures. A similar 
table (table III) was constructed for the load capacity ratio for numerically 
finite. Infinitely short, and Infinitely long bearings. The fluid Is hydrogen. 
This effect remains to be verified experimentally. If It Is correct, this 
demonstrates yet another method of achieving and controlling embedded two-phase 
flow cavities, which do not require a sharp geometrical configuration to 
“trigger" the event. Spontaneous nucleatlon within the metastable field Is 
necessary and sufficient. 

Using the principles of corresponding states, extensions to other single- 
component fluids (perhaps multicomponent fluids) appear reasonable (refs. 13 
and 14) where viscosity becomes 

'.wi-'iH'f)', < 5 > 

and from equation (2) 

». <"- T > • A o( p f- ?) F n < 6 > 

where for single-component fluids 
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for extensions to multicomponent fluids, which are complex (refs. 13 and 14), 
Example : Let the reference fluid (0) be oxygen and the working fluid (a) be 

nitrogen. The bearing temperature and pressure are 100 K and 1.5 MPa. Assum- 
ing constant properties and ©=*■!, solve equation (2) for nitrogen by 
using the thermophysical properties of oxygen. First we need the thermodynamic 
properties at the critical point. 


Fluid 

Tc 

K 

P C > 

MPa 

g/cm 3 

M, 

g/g-mole 

e 

$ 

f 

h 

Fn 

n 2 

°2 

126.3 

154.78 

3.417 

5.082 

0.3105 

.4325 

28.016 
31 .9988 

1 

1 

0.816 

1 .22 

0.74 


f = 126.3/154.78 x 1 = 0.816 ; h - (28.016/31 .9988)(0. 4325/0. 3105) x 1 * 1.22 
F„ = [(31.9988/28.016) 0.816]°- 5 /l .22(2/3) = 0.74 
P0 = P(f/h) * 1.5 (1.22/0.816) = 2.242 MPa 
T 0 » T/f » 100/0.816 = 122.5 K. 

From the thermophysical properties of oxygen (e.g. , ref. 15) the viscosity at 
(P 0 ,To) is 0.936x10-3 g/cm s. The predicted nitrogen viscosity then becomes 

no F Ti 3 0.936x10-3 x 0,74 ■ 0.693x10-3 g/ cm s 

From the thermophysical properties of nitrogen (ref. 15) the viscosity at 
(P,T) Is 0.777x10-3 g/cm s, which represents a 12-percent error. Once the 
bearing geometry and speed are defined, the solution to equation (2) follows 
from equation (5) 


A ( P.T) 

A cyty 


0.74 


to within 12 percent.* Reynolds equation (ref. 11) Is now solved once using 
the bearing parameter for nitrogen at P * 1.5 MPa and T * 100 K and again 
using the bearing parameter predicted for nitrogen from oxygen properties at 
Pq = 2.242 MPa and T = 122.5 K. The reduced circumferential pressure pro- 
files at the center of the bearing, with speed as a parameter, are illustrated 
In figure 5. 


Boilers and Rotors 

A corollary class of nested two-phase flows occurs In turbomachine systems 
as rotating boilers and cooled rotor windings and In a multiplicity of super- 
conducting machinery applications. 


*In reference 11 the solution of the Reynolds equation Is In dimensional 
form and AP = 1 (pressure units). For solutions that are In dimensionless 
form, with fixed P am f, and both Pq and P much larger than P a mb» then 
A/Ag = fFjj/h. 
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A pressure containment vessel, axial flow of a cooling (heating) fluid, 
and rotation of the fluid about an axis parallel to the flow path result In 
centrifugal and Coriolis body forces, which In turn produce pressure and shear 
forces within the fluid. Further, the radial position of the critical pressure 
surface can be adjusted to occur very close to the fluid-vapor Interface. 

For a cylindrical boiler rotating at N-rpm the fluid to be evaporated Is 
forced against the outer wall at a pressure Pwal 1 3 p c* for M > N 
p wall I s greater than P c and vapor moves from the wall to the region where 
P < P c ,; the two-phase region Is nested. Here the Coriolis forces are small, 
but cannot be neglected when, for example, cooling a superconducting motor 
(generator) to prevent the winding from going normal. Again the pressure In 
the winding coolant tube becomes greater than P c for all M > N while the 
pressure at the centerline of rotation Is less than P c ; the system Is weakly 
nested. In this case two-phase flow can occur within the coolant tube and 
become distributed In an unusual axial pattern because of the Coriolis forces. 
In some Instances It may be possible to find a spiral weakly nested region. 

To the dynamlclst the location and extent of the nested regions may become 
Important only after the system Instabilities are encountered. 


SUMMARY OF RESULTS 

Two-phase cavity flow zones have been shown to be embedded within a super- 
critical pressure envelope. For a fixed geometrical configuration these zones 
are established by controlling Inlet stagnation conditions, exit backpressure, 
or combinations. The zones can be "triggered" by discontinuous changes In 
geometry or by significant pressure gradients. 


For high-performance turbomachine seals and possibly bearings, nested two- 
phase zones can be found In the circumferential and axial coordinates. By 
using the Reynolds equation a surface of Incipient two-phase flows was estab- 
lished. The load capacity of the bearing degrades significantly but not before 
an overshoot occurs. It Is speculated that this overshoot "saves" the design. 
With the corresponding-states relations these results can be extended to other 
fluids. These phenomena remain to be demonstrated experimentally. 
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TABLE I. - BEARING PARAMETERS FOR ANALYTICAL SOLUTION 


Geometry: 

Diameter, mm (in) . .. . * . , . , . . , . , 25.4 (1) 

Length, mm (in) . . . . . 25.4 (1) 

Clearance, mm (in) . ............ 0.0254 (0.001) 

Eccentricity ................ 0.97 

Hydrodynamic: 

Inlet pressure, MPa (psia) ................... 13.8 (2000) 

Outlet pressure, MPa (psia) .................. 13.45 (1950) 

Maximum pressure, MPa (psia) ................. 27.16 (3938) 

Minimum pressure, MPa (psia) 0.079 (11.5) 

Inlet temperature, K ( R) ...................... 20 (36) 

Rotational speed, rad/s (rpm) . , * *. . . 7850 (75 000) 


TABLE II. - CAVITY PRESSURE DIFFERENCE 

< P in - P min )/P c AL0N6 CENTERUNE 

[Clearance, 0.076 mm (0.003 in); radius, 38 mm 
(1.5 in); spacer or bearing length, 76 cm 
(3 in); inlet pressure, 1.628 MPa (236 psia); 
exit pressure, 1.524 MPa (221 psia); 
viscosity, 0.127xl0“ 4 Pa-s (0.187xl0" 8 lbf-s).] 


Inlet 

temperature, 

T in 

Dimension- 

less 

eccentricity, 

a 

Rotational speed, 
N, 
rpm 

K 

*R 

10xl0 3 

50xl0 3 

lOOxlO 3 

22.1 

39.8 

0.1 

0.041 

0.043 

0.046 



.5 

.045 

.064 

.087 



.95 

.24 

.646 

— 

28.1 

50.5 

0.1 

0.041 

0.042 

0.044 



.5 

.043 

.056 

.071 



.95 

.18 

— 

* — 


TABLE III. - LOAD CAPACITIES FOR NUMERICALLY FINITE, INFINITELY SHORT, AND INFINITELY LONG BEARINGS 
[Clearance, 0.076 cm (0.003 in); radius, 38 mm (1.5) in); length, 76 cm (3 in); inlet pressure, 
1.628 MPa C236 psia); exit pressure, 1.524 MPa (221 psia); viscosity, Q.127xlG” 4 Pa-s 
(0.187xl0~ 8 lbf-s).j 


Inlet 

temperature, 

T in 

Dimensionless 

eccentricity, 

a 

Rotational 

speed, 

N, 

rpm 

Bearing 

parameter, 

A 

Numerically 

finite 

bearing 

Infinitely 

short 

bearing 

Infinitely 

long 

bearing 

Pinkus/ 

Sternlicht 

K 

KS 

22.1 

39.8 

0.1 

10xl0 3 

5.87 

0.99 

0.71 

4.1 

1.4 




50 

29.3 

4.9 

3.5 

20.6 

7.2 




100 

58.7 

9.9 

7.1 

41.3 

14.3 ' 



.5 

10 

5.87 

6.3 

6.6 

21.2 

2.6 




50 

29.3 

31.3 

33.1 

106 

13.2 




100 

58.7 

62.7 

66.1 

212 

26.4 



.95 

10 

5.87 

53 

863 

87 

20 




50 

29.3 

265 

4318 

433 

100 




100 

58.7 


— 

— 

200 

28.1 

50.5 

0.1 

lOxlO 3 

3.9 

0.66 

0.47 

2.7 





50 

19.5 

3.3 

2.3 

13.7 





100 

39.1 

6.6 

4.7 

27.5 

— 



.5 

10 

3.9 

4.2 

4.4 

14.1 

HM „n . — m. — 




50 

19.5 

20.9 

22 

70.6 

. 




100 

39.1 

41.7 

44 

141 

— 



.95 

10 

3.9 

35.4 

574 

57.7 

— 
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Figure L- Temperature-entropy schematic. 
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(b) Applied backpressure. 

Figure Z - Pressure profile for fully eccentric, nonrotating, three-step 
SSME seal configuration. 
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Figure 4.- Pressure profile for nitrogen flow through high-L/D tube 
with Borda inlet. (Dimensions are in centimeters (inches).) 
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Figure 5. - Reduced midbearing circumferential pressure profiles 
based on corresponding states for fluids nitrogen and oxygen. 
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